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Abstract 


Using data collected during spring and summer seasons in the north polar Atlantic we examined 
the variability of the spectral absorption, a(k), and backscattering, b\,(k), coefficients of surface 
waters and its relation to phytoplankton pigment concentration and composition. For a given 
chlorophyll a concentration ( TChla ), the concentrations of photosynthetic carotenoids (PSC), 
photoprotective carotenoids (PPC), and total accessory pigments (AP) were consistently lower in 
spring than in summer. The chlorophyll-specific absorption coefficients of phytoplankton and 
total particulate matter were also lower in spring, which can be partly attributed to lower 
proportions of PPC, PSC, and AP in spring. The spring values of the green-to-blue band ratio of 
the absorption coefficient were higher than the summer ratios. The blue-to-green ratios of 
backscattering coefficient were also higher in spring. The higher bh values and lower blue-to- 
green hb ratios in summer were likely associated with higher concentrations of detrital particles 
in summer compared to spring. Because the product of the green-to-blue absorption ratio and the 
blue-to-green backscattering ratio is a proxy for the blue-to-green ratio of remote-sensing 
reflectance, we conclude that the performance of ocean color band-ratio algorithms for 
estimating pigments in the north polar Atlantic is significantly affected by seasonal shifts in the 
relationships between absorption and TChla as well as between backscattering and TChla. 
Intriguingly, however, fairly good estimate of the particulate beam attenuation coefficient at 660 
nm (potential measure of total particulate matter or particulate organic carbon concentration) can 
be obtained by applying a single blue-to-green band-ratio algorithm for both spring and summer 
seasons. 
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1. Introduction 


Empirical approaches, in which in situ oceanographic measurements are used to relate band 
ratios of ocean spectral reflectance to surface water properties such as chlorophyll a 
concentration ( Chla ), are an important component of routine processing of ocean color satellite 
imagery [e.g., O'Reilly et al., 1998], The well-known example is the OC4 (Ocean Chlorophyll 4) 
empirical algorithm that is used for processing of global satellite data collected with SeaWiFS 
(Sea-viewing Wide Field-of-view Sensor). This algorithm estimates chlorophyll a concentration 
from a fourth order polynomial by using the maximum band ratio determined as the greatest of 
the i? rs (443)/f? rs (555), f? rs (490)/f? rs (555), or R rs (5 10)/R rs (555) values [ O'Reilly et al, 2000]. R n (X) 
is the spectral remote-sensing reflectance defined as the ratio of the spectral upwelling water- 
leaving radiance just above the sea surface to spectral downwelling plane irradiance incident on 
the sea surface, where X is a wavelength of light in vacuum. 

A basic connection between 2-band algorithms and inherent optical properties (IOPs) of 
seawater comes from an approximate relation [ Gordon et al., 1988]: 

RM ~ b h (X)/a(X) (1) 

where a(X) and b\iX) are the spectral absorption and backscattering coefficients of seawater, 
respectively. From Equation (1), the reflectance band ratio is approximately equal to the product 
of backscattering ratio and absorption ratio: 

R rs (h)/ROXi) ~ [b h (Xi)/b h (X 2 )] [a(X 2 )la(X x )\ (2) 

which indicates that changes in R rs (Xi)/R rs (X 2 ) are driven primarily by the variability in spectral 
ratios of bb(Xi)/bb(X 2 ) and a(X-i)la{X\). Therefore, concurrent measurements of b\XX) and a(X) 
offer the opportunity to explain how these IOPs affect the general trends and the spread of data 
in the relationship of chlorophyll a concentration versus /? rs ( A i )/'/? rs ( /_ 2 ) as well as the temporal 
and geographical variability in this relationship. The a(X) and b\,(X) coefficients are usually 
partitioned into a few additive components associated with broadly defined categories of 
seawater constituents [e.g., Mobley, 1994]: 

a(X) = a w ( X) + a p (X) + <7cdom(A.) (3) 
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a p (X) = a v lh (X) + a d (X) 
b b (X) = b bw (X) + b bp (X) 

where the subscript w denotes pure seawater, p the total particulate assemblage suspended in 
water, ph the phytoplankton suspended in water, d the detrital or non-phytoplankton particles 
(this component in fact includes non-living organic and mineral particles, and heterotrophic 
organisms), and CDOM the colored dissolved organic matter. The partitioning of IOPs into 
these components offers an added value for the interpretation of the effects of b b (X\)/b b (X 2 ) and 
a(X 2 )/a(X\) on the relationship between Chla and R rs (Xi)/R rs ( X 2 ). 

A limited amount of appropriate in situ data has held back a development of regional bio- 
optical and ocean color algorithms in the high northern latitude waters, where relatively few 
experiments with a broad suite of bio-optical measurements have been made so far [Mitchell, 
1992; Cota, 2001; Sathyendranath et al., 2001]. Recently, we described bio-optical relationships 
derived from data collected on R/V Oceania during the summer season (June - August) in the 
eastern part of north polar Atlantic between Norway and Spitsbergen [Stramska et al., 2003]. 

The main goal of the present study is to compare these previous results with similar relationships 
established from more recent data collected on R/V Polarstern during the spring season (April- 
May). These latest data were collected on a ship’s transect starting near Spitsbergen and 
extending west all the way to the vicinity of Greenland. The eastern part of the Polarstern 
transect overlapped, in part, with the study area investigated earlier during the Oceania summer 
cruises. We discuss the variability of the spectral absorption and backscattering coefficients of 
surface waters and the relation between these IOPs and surface chlorophyll a concentration in the 
north polar Atlantic. This analysis provides insights into seasonal and spatial variability of bio- 
optical properties that are essential to ocean color algorithms for estimating Chla from the blue- 
to-green reflectance ratio within the investigated region. 
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2. Data and methods 


Data were collected in the north polar Atlantic during spring and summer seasons. The 
spring data were collected in April-May 2003 on R/V Polarstern operated by Alfred Wegener 
Institute for Polar and Marine Research (AWI) in Gennany. The summer data were collected in 
June- August of 1999 and 2000 during two cruises on R/V Oceania operated by Institute of 
Oceanology, Polish Academy of Sciences (IOPAS). Whereas the Oceania cruises covered the 
eastern part of the north polar Atlantic (0°-20°E and 70°-80°N), the Polarstern transect extended 
between 10°E and 16°W, mostly along 75°N. The maps showing the location of stations can be 
found elsewhere [ Stramska et al., 2003; Stramska and Stramski, 2005]. The data from Polarstern 
were divided into two subsets. First, a subset #1 includes data from the eastern side of the 
transect (between Spitsbergen and 0°E, mostly at 75°N), and second, a subset #2 includes data 
from the western side of the transect (between 4°W and Greenland, mostly at 75°N). This 
division is justified by HPLC analysis of phytoplankton pigments, which suggests significant 
differences in phytoplankton composition on both sides of the transect, as discussed below. In 
addition, the spring Polarstern data of subset #1 were collected in the geographic area 
overlapping with some summer stations occupied on the Oceania cruises. This facilitates 
unambiguous comparison of data for seasonal variability. However, comparisons of all data 
from Polarstern and Oceania cruises reflect not only seasonal but also spatial variability within 
the entire study region. 

In this study we consider data from stations where we observed no evidence of significant 
influence of terrestrial material on the optical properties of seawater. Because our interest is on 
bio-optical relationships relevant to remote sensing of ocean color, we use data representing 
surface waters only. During the Polarstern cruise in spring, surface waters in the north polar 
Atlantic were well-mixed or weakly stratified. Therefore we assume that data collected between 
the surface and a depth of 30 m characterize the surface water properties and can be used for our 
purposes. In contrast, during the Oceania cruises in summer, surface waters were stratified, so 
we use data collected only between the surface and 10 m. In this paper, we do not discuss in situ 
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radiometric measurements. Although such data were collected on Oceania [ Stramska et al., 
2003], the cruise schedule and ice conditions allowed only for a few underwater radiometric 
measurements from Polarstem, which is not sufficient to make a meaningful comparison of the 
radiometric data between the cruises. 

Methods of our bio-optical measurements that were carried out on R/V Oceania and R/V 
Polarstem are described in Stramska et al. [2003] and Stramska and Stramski [2005]. Below is a 
description of methods of direct relevance to this paper. 

2.1. Pigments 

For the analysis of phytoplankton pigments, suspended particles were collected by fdtration 
of water samples onto Whatman glass-fiber filters (GF/F) under low vacuum. Three methods 
were used to measure the concentration of pigments. The spectrophotometric method (UV4-100 
spectrophotometer, Unicam, Ltd) was used to estimate chlorophyll a concentration ( Chla ) on 
board the ship during the Oceania cruises in 1999 and in 2000 (77 samples). In this method the 
optical density (absorbance) of pigment extract in ethanol was measured at 665 nm. After 
correction for background signal in the near infrared (750 nm), the absorbance was converted to 
chlorophyll a concentration, Chla, using the chlorophyll ^/-specific absorption coefficient in 96% 
ethanol [see Stramska et al., 2003 for details]. In addition, pigment samples (3 1 samples) on the 
Oceania cruise in 1999 were stored in liquid nitrogen and analyzed with fluorometric (Chla and 
phaeopigments concentrations) and high-performance liquid chromatography (HPLC) methods 
after the cruise [ Bidigare and Trees, 2000; Trees et al., 2000a], Fluorometric and HPLC 
methods were also used to analyze pigment samples collected on the Polarstem cruise. Here we 
report on data from the analysis of 76 pigment samples collected in the top 30 m of water on 
Polarstem. The HPLC analysis provided an estimate of total Chla ( TChla ) concentration that 
was calculated as a sum of chlorophyll a and derivatives (chlorophyllide a, chlorophyll a 
allomers and epimers). HPLC also provided concentration of various accessory pigments. 
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The comparison of chlorophyll a concentrations obtained with different methods is shown in 
Figure 1 . The correlation of various estimates is quite good. In particular, a good agreement is 
observed between the spectrophotometric Chla and HPLC-derived TChla for samples from the 
Oceania cruise in 1999. Because of this agreement and because we did not carry out HPLC 
analysis on Oceania in 2000, we will use the spectrophotometric Chla estimates as the 
chlorophyll a concentration for both Oceania cruises. When we discuss the HPLC data from 
Oceania, it will be explicitly indicated. For the Polarstem cruise we will use the HPLC-derived 
TChla estimates for the chlorophyll a concentration. 

For more detailed comparisons of phytoplankton pigments from HPLC analysis between our 
data sets, we define three major categories of accessory pigments. The first category, the 
photosynthetic carotenoids (PSC), includes peridinin (per), fucoxanthin (fuco ), 
19’hexanoyloxyfucoxanthin (hex), and 19’butaonoyloxyfucoxanthin (but). The second category, 
the photoprotective carotenoids (PPC), consists of diadinoxanthin (dda), alloxanthin (alio), 
diatoxanthin (diato), zeaxanthin (zea), violaxanthin (viol), and a- and [1-carotcnes. The third 
category, referred to as the total accessory pigments (AP), is defined as the sum of PSC, PPC, chi 
b, and chi c. 

HPLC analysis of phytoplankton pigments provides data that can be used to characterize 
broad algal classes present in a water body. One of the problems in the use of HPLC data for 
describing phytoplankton community structure is the fact that algal groups do not necessarily 
have unique pigments to allow unambiguous taxonomic identification. Various approaches have 
been tested over the years to overcome this problem. One of the most successful methods is 
based on factor analysis developed into a computer program CHEMTAX [Mackey et ah, 1996]. 
This program was originally used to estimate algal class abundance for samples collected 
between Australia and Antarctica, where the predominant phytoplankton cells in terms of 
contribution to Chi a are relatively large [ Wright et al., 1996] and in the western equatorial 
Pacific where picoplankton dominate [Mackey, et al., 1998]. Comparisons with traditional cell 
counting techniques supported the results obtained with CHEMTAX for these locations. Other 
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successful applications of CHEMTAX include studies of phytoplankton community structure in 
the Canadian Arctic [ Vidus si et a/., 2004] and near Faroe Islands [Riegman and Kraay, 2001]. 

Using our HPLC data we applied CHEMTAX to estimate the contributions of various algal 
classes to the TChla in our study region. As input to the code, CHEMTAX uses the accessory 
pigment/TC/z/a ratios for a number of algal classes. We did not have available information on 
accessory pigment/TC/z/o ratios for phytoplankton classes present in our study region. Therefore 
we used the pigment fractions shown in Table 1, which are similar to those used recently to 
examine phytoplankton community structure in the Faroe-Shetland Channel located to the south 
of our study region [Riegman and Kraay, 2001]. These fractions are in general agreement with 
other data found in the literature [Jeffrey and Wright, 1997]. The fractions given in Table 1 were 
iteratively modified by CHEMTAX to minimize the sum of squared differences between the 
observed and calculated pigment concentrations. Because pigment ratios are known to vary with 
environmental growth conditions (light, nutrients, temperature), we carried out separate 
calculations for data collected on Oceania and Polarstern cruises to reveal possible seasonal and 
regional differences. In addition, because our initial CHEMTAX runs suggested different 
phytoplankton composition at the eastern (data subset #1) and western (data subset #2) parts of 
the Polarstern transect, independent CHEMTAX runs were made for each of these Polarstern 
data subsets. 

The HPLC data were also used to estimate proportions of phytoplankton present in broad size 
groups using the method developed by Vidussi et al. [2001], According to this method, if the 
total diagnostic pigments (DP) are defined as the sum of Chi h, zea, alio, but, hex,fuco, and per, 
then the relative contribution to TChla by picoplankton (< 2 pm), nanoplankton (2-20 pm) and 
microplankton (> 20 pm) can be estimated from the following relationships: 

Picoplankton (%) = (zea + Chi b) ’ 100/DP 
Nanoplankton (%) = (alio + but + hex) ' 100/DP 
Microplankton (%) = (fuco +per ) ' 100/DP 
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where all pigment concentrations are in mg m" . Derived in this way phytoplankton size 
fractions present some limitations [Vidussi et al, 2001]. These limitations result from the facts 
that certain pigments are present in more than one taxonomic group and that the same 
phytoplankton taxa can be generally present in different size classes. Nevertheless, to first 
approximation, we can use this tool to test whether there are significant differences in size 
structure of phytoplankton populations within the north polar Atlantic. 

2.2. Inherent optical properties 

In situ vertical profiles of IOPs and physical properties of seawater were measured with a 
multisensor datalogger system that consisted of the SeaBird Sealogger 25 and SeaBird 
temperature, conductivity, and pressure sensors, Hydroscat-6 backscattering sensor (442, 470, 
555, 589, 620, and 671 mn, HobiLabs, Inc.), and two c-Star beam transmissometers (488 and 660 
nm; WetLabs, Inc). We processed the Hydroscat-6 data with the Hydrosoft software (version 2.6 
of December 2002, HobiLabs, Inc.), which includes corrections suggested by Boss and Pegau 
[2001]. Because the c-Star transmissometer data were occasionally unavailable, we consistently 
used a standard correction built-in the Hydrosoft software to account for the light attenuation 
effects. We verified that the use of correction based on actual beam attenuation measurements 
instead of the standard correction generally yielded the b b (k) values higher by a few percent but 
the blue-to-green ratios of b b were nearly the same (to within 2% in the extreme cases). 

Therefore, the choice of attenuation correction method for determining b b does not affect main 
conclusions from our analysis. Our final b b (X) values at Hydroscat-6 wavebands were used to fit 
a power function, b b (X) ~X ~ y , which describes the spectral shape of b b (X) spectrum. The b b (X) 
values at wavelengths corresponding to nominal center wavelengths of SeaWiFS bands were 
then obtained from this function. 

During processing of optical data, all profiles were first carefully inspected for quality, for 
example for the presence of possible noise in the near-surface data. Doubtful or noisy data were 
removed from the analysis. Typically we did not consider data acquired within the top 2 or 3 
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meters of the ocean. The remaining data were averaged into 2-m bins to provide final depth 
profiles. Surface IOPs were approximated as the averages for the uppermost bin. The 
relationships for Chla versus bh(k) and Chla vs. c p (660) were obtained by matching Chla 
estimates from discrete water samples and the IOPs measured at corresponding depths. The time 
difference between the collection of discrete water samples and acquisition of in situ optical data 
was usually less than an hour. 

The spectral absorption coefficient of particulate matter, a p (k), was derived from filter-pad 
measurements using a dual-beam UV4- 100 spectrophotometer (Unicam, Ltd). First, an 
appropriate volume of seawater (0.5 - 2 L depending on particle concentration in water) was 
filtered onto a 25-mm glass-fiber filter (GF/F). Next, the absorption spectrum of particles, a p (k), 
retained on a filter was determined using the transmittance-reflectance (T-R) spectrophotometric 
method [Tassan and Ferrari, 1995 a, b]. Transmittance and reflectance were measured between 
380 and 750 nm with a 1 nm interval. Correction for pathlength amplification factor (P) was 
made using the following expression [Kaczmarek et ah, 2003]: 

OD s = 0.592 ODf 2 + 0.4 OD f 

where OD s is the optical density (absorbance) of particle suspension with no p effect and ODf is 
the measured optical density of particles on the filter. 

The absorption coefficient of non-phytoplankton (detritus) particles, ad(T), was determined 
with another run of transmittance and reflectance measurements after bleaching the filter with 
sodium hypochlorite [Ferrari and Tassan, 1999]. Because it was recently shown that oceanic 
particles generally do not absorb significantly in the near infrared spectral region [Babin and 
Stramski, 2002], the final estimates of a p (k) and adk) were derived after subtracting the 
measured values of a p and a& averaged between 740 and 750 nm. Finally, the phytoplankton 
absorption coefficient, a p h(k), was calculated as a difference between a p (k) and ad'k)- During the 
Oceania cruises the absorption measurements were made on freshly collected samples on board 
the ship. The a p samples from the Polarstern cruise were stored in liquid nitrogen and analyzed 
in the laboratory after the cruise. 
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The absorption data from the Oceania cruises presented in our previous paper [Stramska et 
al., 2003] were calculated from filter-pad measurements performed in the transmittance (T) 
configuration [ e.g ., Mitchell et al., 2000]. However, because the reflectance (R) measurements 
were also made on Oceania, we now recalculated the absorption spectra using the T-R algorithm. 
In the blue-green spectral region the differences between the ’old’ and 'new' a v (X) estimates from 
Oceania cruises are generally within a few percent. We note that the T-R method may be 
superior to T-method. The T-R method accounts (albeit still not perfectly) for wavelength- 
dependent scattering error and the pathlength amplification factor applied in the T-R method is 
less variable than in the T-method [Ferrari and Tassan, 1999; Tassan and Ferrari, 2002]. The 
absorption data from the Oceania and Polarstem cruises discussed in this paper were detennined 
in the same way with the T-R method. 

3. Results and discussion 
3.1. Pigments 

Surface chi a concentrations varied from 0.2 to 0.8 mg m~ 3 and from 0.2 to 10 mg m' 3 in the 
eastern and western parts of the Polarstern transect respectively. High chi a concentrations in the 
west indicate the presence of a well-developed phytoplankton bloom, while comparatively low 
chi a concentrations in the east suggest a pre-bloom or early bloom phase. Fairly high chi a 

•3 

concentrations observed during Oceania cruises (0.2 - 3 mg m" ) suggest that phytoplankton 
growth in the eastern part of the region has been supported throughout the summer season, most 
likely due to periodical replenishment of nutrients by mixing events caused by episodic storms. 

CHEMTAX pigment analysis revealed some differences in the composition of phytoplankton 
community between our data sets, that is between the summer data set from Oceania and two 
spring subsets (eastern and western) from Polarstern. Note that as explained in Mackey et al. 
[1998] the phytoplankton classes derived from CHEMTAX do not correspond exactly to 
conventional taxa but are representative of typical pigment composition. 
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In our data sets, diatoms and prymnesiophyceae generally dominated the phytoplankton 
biomass expressed in terms of TChla. At nearly all stations, diatoms and prymnesiophyceae 
combined together contributed more than 50% of TChla (Figure 2). The dominance of diatoms 
was most pronounced in the western part of the Polarstern transect, where we observed an active 
development of spring phytoplankton bloom. Stations located in the eastern part of the 
Polarstern transect, which likely represent a pre -bloom situation, were dominated by 
prymnesiophyceae (~ 30% of TChla). Other phytoplankton groups with significant contribution 
to TChla in the eastern part of Polarstern transect were chrysophyceae and prasinophyceae (each 
~ 20% of TChla). The stations visited in summer by Oceania were dominated either by diatoms 
or by prymnesiophyceae (on average each group represented somewhat less than 30% of TChla). 
Cryptophyceae and chrysophyceae each contributed on average a little more than 10% of TChla 
during Oceania and Polarstern cruises. The summer stations were characterized by somewhat 
more abundant populations of dinophyceae and cyanophyceae and less abundant population of 
prasinophyceae than the spring stations, but the contribution to TChla of each of these groups 
was on average low (well below 10%). The chlorphyceae were present in a small percentage of 
the total phytoplankton community. Examples of the final pigment ratios estimated with 
CHEMTAX are shown in Figure 3. The high ratios for fuco and hex for diatoms and 
prymnesiophyceae in Oceania data set are particularly noticeable. 

In Figure 4 the concentrations of three pigment groups (PPC, PSC, and AP) are plotted as a 
function of TChla. Each pigment group is well correlated with TChla. The most important 
observation is that the concentrations of PPC, PSC, and AP for a given TChla were consistently 
higher in summer (Oceania data set) than in spring (both for the eastern and western data from 
Polarstern). For example, the concentrations of AP for a given TChla were, on average, about 
two times higher in summer than in spring. The parameter values for the best fit of the linear 
regression on log-transfonned data for total accessory pigments versus TChla are compared in 
Figure 4c with similar global relationship established by Trees et al. [2000b]. In contrast to our 
study, Trees et al. included phaeopigments in the accessory pigments. However, their results of 
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statistical analysis were not significantly affected by phaeopigments because of low 
concentrations of these pigments in their data set. Our data are within the range reported by 
Trees et al. However, most of our data points from the summer season are above the global 
relationship of Trees et al. and most of our spring data are below that relationship. 

Figure 5 shows proportions of pico-, nano-, and microphytoplankton derived from the 
method of Vidussi et al. [2001], Picoplankton was always the least important of the size classes 
in terms of contribution to TChla. The lowest and the highest relative contributions to total 
phytoplankton biomass from microphytoplankton were detected in the eastern and the western 
parts of the Polarstern transect, respectively. Generally, the differences in phytoplankton size 
structure seem to be more pronounced between the eastern and western sections of the Polarstern 
transect, than between the averages for the entire Polarstern and Oceania data sets (the averages 
for combined Polarstern data are not shown in Figure 5). 

3.2. Phytoplankton and total particulate absorption 

Phytoplankton were the dominant component of particulate absorption for both the Polarstern 
and Oceania data sets as shown by the relationships between a P h(A.) and a v (X) for the blue 
wavebands, 442 and 490 nm in Figure 6. In spring, phytoplankton contributed, on average, 76% 
and 79 % to a p at 442 and 490, respectively. In summer this contribution was about 71% and 
72%, respectively. This indicates that the percent contribution of detritus to absorption in the 
blue was generally somewhat higher in summer than in spring. In the green spectral band 
centered at 555 nm, phytoplankton contributed, on average, about 64% to the total a p in both the 
spring and summer seasons (not shown). 

Figure 7 shows representative examples of <7 p h(a) nonnalized to chlorophyll a concentration, 
i.e., the chlorophyll-specific absorption coefficient of phytoplankton, a p h(k). This figure 
illustrates a trend of decrease of fl* p h(A.) with an increase of chi a concentration. Similar trend 
was observed in larger data sets that combined observations from various regions [Bricaud et al., 
1995; Babin et al., 2003] and it was attributed to increasing pigment packaging, changes in 
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species composition, and changes in proportions of accessory pigments from oligotrophic to 
eutrophic waters. Importantly, Figure 7 also shows a significant difference between the 
relatively high a p h(A.) observed in summer on the Oceania cruises and low values in spring on 
the Polarstern cruise. Note that both the eastern and the western parts of the Polarstem transect 

m sH ... 

were characterized by lower a p h(A.) than the summer values from the Oceania cruises. This 
indicates a seasonal variation in a p h(X). This variation can be, to a large degree, attributed to 
higher proportions of PPC, PSC, and AP in summer than in spring populations of phytoplankton 
(see Figure 4). This trend of seasonal change in the ratio of accessory pigments to TChla is 
consistent with similar seasonal patterns in pigment composition described by others [Babin et 
al., 2003]. 

The difference between the spring and summer data of a P h*(A.) can be also partly attributed to 
the pigment package effect. The package effect reduces the a v \* values and flattens the a v \*0-) 
spectrum. This effect increases with an increase of average cell size and intracellular 
concentration of pigments within phytoplankton cells [Morel and Bricaud, 1981]. The package 
effect is often estimated by comparison of a p h*(675) with the pigment absorption of chlorophyll 
a in solution because accessory pigments have relatively small influence on a p h* in the red 
absorption band of chlorophyll a. Because a p h*(675) was, on average, 0.17 and 0.27 m (mg chi 
a)~' for Polarstem and Oceania data sets, the package effect appears to be stronger for the spring 
phytoplankton populations than for the summer populations. Based on the differences in 
phytoplankton size structure shown in Figure 5, one could also expect a regional differentiation 
in the package effect between the eastern and western Polarstern data. In the western part of the 
Polarstern transect the package effect could have been higher than in the eastern part due to a 
higher contribution of large microphytoplankton cells to the total phytoplankton biomass as 
indicated in Figure 5. However, no clear differences in the magnitude of a p h* are observed 
between the two Polarstern data subsets in Figure 7. We note that variations in the concentration 
of chlorophyll b and phaeopigments also contribute to a p h*(675), which may confound the 
interpretation of the package effect. 
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Because both the eastern and western subsets of Polarstem data show significantly lower 
a ph(A-) than the Oceania data, from now on we will discuss and compare the differences in the 
absorption between the two main seasonal data sets, that is the summer data set (from Oceania) 
and the spring data set (i.e., combined eastern and western Polarstem data). We will, however, 
keep on distinguishing the eastern and western Polarstem data points on most graphs. The 
differences in absorption between the summer and spring data sets are reflected in Figure 8, 
which shows the relationships for a p h and a p at 442 and 555 mn vs. chi a concentration. The a p h 
values for a given chi a concentration are consistently lower in spring than summer (Figure 8a). 
Similarly, the a p values are generally lower in spring than in summer (Figure 8b). Therefore, 
similarly to a p h*, the chlorophyll ^-specific absorption coefficient for total particulate matter, 
a p *, is lower in spring than in summer. Figure 8 also compares our absorption data with 
relationships obtained by Bricaud et a/. [1995; 1998] for a large data set from various regions. 
Whereas our summer data tend to be higher than the values predicted with the Bricaud et al. 
regression lines, our spring data tend to be lower than these lines. 

Differentiation between our spring and summer data sets from the north polar Atlantic is also 
obvious if we compare the green-to-blue spectral ratios of the absorption coefficient vs. pigment 
concentration (Figures 9, 10, and 11). If we consider a sum of the absorption contributions of 
pure water and phytoplankton, a w+ph , there is a good correlation between the a w+p h band ratio and 
chi a concentration if each of our data sets is considered separately. However, for the same chi a 
concentration, the spring data show consistently higher values of <7 W + P h band ratio than the 
summer data (Figure 9). Similar behavior of the a w + P h band ratio is observed when plotted vs. 
total pigment concentration (i.e., vs. TChla + AP) although in this case the seasonal differences 
are somewhat smaller (Figure 10). When we consider a sum of absorption contributions of pure 
water and total particulate matter, a w+p , the seasonal differences in the a w+p band ratio at any 
given chi a concentration are even more pronounced than in the case of the n w + P h band ratio 
(Figure 1 1). The seasonal difference in the ratio of TChla to accessory pigments (Figure 4) 
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appears to be a significant cause for the seasonal differences in the data shown in Figures 9 and 

11 . 

The behavior of the green-to-blue absorption ratio has important implications for ocean color 
algorithms that are based on the blue-to-green reflectance ratio (see Equation 2). It is important 
to realize that our absorption band ratios in Figure 1 1 do not fully represent the absorption term 
on the right-hand side of Equation 2. This is obviously because we have omitted the contribution 
of colored dissolved organic matter, <2cdom(^), to the total absorption coefficient, a(X). 
Unfortunately, ocdom(A-) was not measured during our experiments. However, the contribution 
of acDOM(^) to a(k) was suggested to be quite small for our summer data set \Stramska et al., 
2003]. For open waters within the investigated region, we expect that «cdom(^) would not play a 
significant role in spring either because the estimates of DOM concentration during the 
Polarstem cruise were relatively low ( Schwarz et al., submitted). Therefore, it is reasonable to 
assume that the green-to-blue ratio of the total absorption coefficient would behave qualitatively 
in a similar manner as the green-to-blue ratio of the absorption coefficient by pure seawater and 
particles displayed in Figure 11. 

3.3. Backscattering 

The backscattering coefficients b b (k) at 442 and 555 nm for any given chi a concentration are 
generally lower for the Polarstern than for the Oceania data set, although the scatter of data 
points is large, especially for the Oceania data set (Figure 12a). In spring the contribution of 
particles to total h|,(555) (i.e., the sum of particulate and molecular backscattering at 555 nm) 
was usually between 40 and 50%, and exceeded 60% only at a few stations. In summer this 
contribution at most stations was between 70 and 90%. The spectral dependency of hb(k) was 
also different for the two data sets (Figure 12b). Assuming that b b (X) has the wavelength 
dependence defined by: 

b b (X) = [h bw (555) + h b p(555)](A,/555)" T (4a) 

the slope y for the spring data is best described by the following equation: 
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(4b) 


y = 0.026 6 b (555)' 0 ' 750 

and for the summer data by: 

y = 0.067 6 b (555)" 0 ' 572 (4c) 

The changes in y result primarily from variability in relative contributions of water molecules 
and particles to the total backscattering [Morel and Gentili, 1991]. Pure water backscattering 
coefficient, 6 bw (k), is characterized by strong wavelength dependency, h bw ( X) ~ X' ' [Morel, 
1974]. The wavelength dependency of particulate backscattering coefficient, 6 bp ( X), for natural 
populations of marine particles is usually less pronounced, and 6 bp (k) ~ /." 1 is sometimes 
assumed in models [Carder et al., 1999]. In general, however, 6 bp (k) depends on the variability 
in concentration, size distribution, and composition of particulate assemblages, and these sources 
of backscattering variability are not yet well understood [Stramski et al, 2004], Note that in 
Figure 12b the steepest slopes y close the value of 4 or even slightly exceeding 4 were observed 
in clear waters with low 6 b (555) (corresponding to relatively low chi a concentration) in the 
spring season. In these cases the water molecules dominated the total backscattering coefficient. 

According to Equation (2), the blue-to-green backscattering ratio is one of the important 
factors driving the reflectance band ratio that is used in chlorophyll algorithms. Figure 13 shows 
6 b (442)/h b (555) and 6 b (490)/6 b (555) plotted as a function of chi a concentration for the spring 
and summer data sets. For each of the data sets, chi a concentration changes nearly 30-fold but 
the range of h b ratios is relatively small, that is a factor of 1.6 and 1.3 for 6 b (442)/6 b (555) and 
6 b (490)/6 b (555), respectively. Most importantly, however, these backscattering ratios are 
significantly higher for the spring data than the summer data, which is consistent with generally 
higher spectral slopes in the spring ( see Figure 12b). The results in Figure 13 have important 
implications for the seasonal differentiation of the band ratio algorithms for estimating chi a 
concentration. The differences in the h b band ratios shown in Figure 13 will reinforce the effect 
of absorption because the green-to-blue absorption ratio (Figure 1 1) was also higher in spring 
than in summer. 
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Without detailed data of particle size distribution and composition, the interpretation of 
the observed differences in backscattering would be mostly speculative. It was, however, 
hypothesized that non-living particles (mostly from the small size range < 1 pm) can dominate 
the particulate backscattering in most open ocean situations in the absence of intense 
phytoplankton blooms [Morel and Ahn, 1991; Stramski and Kiefer, 1991; Stramski et a/., 2004]. 
Our absorption data suggest that non-living particles were more abundant in summer than in the 
spring season because we observed higher total particulate and higher detrital absorption 
coefficients as well as higher POC concentrations for a given chi a concentration [not shown 
here] during summer. Therefore, higher bb values and lower blue-to-green 6 b ratios observed in 
summer can perhaps be associated with higher concentrations of detrital particles observed in 
summer compared to spring. The backscattering ratio of particles defined as 6 bp (k)/6 p (k) where 
b p (k) is c p (X) - a v (X), is sometimes used as a source of general infonnation about particle sizes 
and/or composition [see Stramski et al., 2004 and references therein]. We found that at A = 488 
mn, the 6 bp /6 p values in summer were either lower or similar to the values measured in spring 
(Figure 14a). In contrast, at X = 660 mn, the 6 bp /6 p values in summer were either higher or 
similar to the values measured in spring (Figure 14b). We show the results for these two 
wavelengths because the beam attenuation c p (k) was measured at these two bands only. One 
conclusion from Figure 14 is that the 6 bp /6 p values are relatively small (usually less than 0.01), 
which suggests that particulate assemblages in the investigated region are dominated by organic 
particles with low refractive index. Another important conclusion from Figure 14 is that the 
particulate assemblages in spring were characterized by much stronger wavelength dependency 
of the backscattering ratio than in the summer season. In spring, A hp (48 8 )/A p (488 ) is generally 
greater than 0.004 - 0.005 but 6 bp (660)/6 p (660) is always smaller than 0.004 and often times less 
than 0.002. For the summer data set, the spectral dependency of the bh p /b p ratio is not so obvious 
or well pronounced. These differences in the spectral behavior of 6 bp /6 p suggest that particles in 
summer could have been, on average, smaller than particles in spring [Morel and Bricaud, 1986; 
Stramski et al. 2001], In addition, the overall particle concentration in summer was probably, on 
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average, higher than in spring, as suggested by higher beam attenuation in summer ( see section 
3.4 below). The possible seasonal differences in particle concentration and size distribution 
would support generally higher backscattering per unit chlorophyll a concentration in summer 
than in spring shown in Figure 12. 

3.4. Particulate beam attenuation and single scattering albedo 

The relationship between the particulate beam attenuation coefficient at 660 mn, c p (660), and 
chlorophyll a concentration shows a shift between the spring and summer data (Figure 15). For 
any chi a concentration, c p (660) is generally higher during the summer season than in spring. 

This shift is qualitatively similar to that discussed above with regard to the particulate absorption 
and backscattering coefficients. The seasonal shifts in the IOP vs. chi a concentration 
relationships appear to be related to increased roles of accessory pigments and/or detrital 
particles in summer. The variability in particulate absorption, especially by accessory pigments, 
is expected to have minor influence on the particulate beam attenuation at 660 mn. Therefore the 
observed seasonal shift in c p (660) vs. chi a concentration is most likely due to an increased 
detrital scattering in summer. 

Given significant seasonal differences in several bio-optical relationships discussed above, it 
is remarkable that a range of single scattering albedo values, co op (k), for suspended particulate 
matter is about the same for the spring and summer seasons (Figure 16). The co op (A) values 
calculated as 6 p (k)/c p (k) where b v (X) = c p (a) - a p (A,), characterize the particles in terms of the 
relative roles of scattering and absorption in the beam attenuation process. Figure 16 indicates 
that scattering dominates the beam attenuation both in the blue band at 488 mn and in the red 
band at 660 mn. Figure 16 also suggests that m„ p (a) at any given phytoplankton biomass (chi a 
concentration) shows no clear systematic differences between the spring and summer seasons. 

One of the most striking results is that the green-to-blue absorption ratio follows one 
consistent pattern as a function of c p (660) regardless of whether the data were collected in spring 
or summer (Figure 17). This is not the case for the blue-to-green backscattering ratio, which 
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shows a seasonal shift if plotted vs. c p (660) (Figure 18). This shift is qualitatively similar to that 
observed for the IOP band ratios plotted as a function of Chi a (see Figures 1 1 and 13). 

3.5. Implications for ocean color band-ratio algorithms 

We now illustrate variations of chi a concentration, total pigments (= TChIa+ AP), c p (660), 
and a p h(442) as a function of the product of the backscattering ratio hb(442)/hb(555) and the 
absorption ratio a w + P (555)/a w + P (442) (Figure 19). Similar illustration but for the 490 nm and 555 
nm bands is shown in Figure 20. Because the product of the IOP band ratios is a proxy for the 
blue-to-green ratios of remote-sensing reflectance R vs (see Equation 2), these figures show 
potential variability in the empirical band-ratio algorithms for estimating pigments, 
phytoplankton absorption, and particulate beam attenuation in the north polar Atlantic. We 
recall a limitation of these results due to the fact that the presented absorption ratios do not 
include a contribution from colored dissolved organic matter. Nevertheless, we expect that the 
main effect of CDOM would be to shift slightly the presented patterns of data points to the left 
without much effect on the differentiation between the spring and summer data sets. 

The results in Figures 19 and 20 are straightforward consequences of the relationships 
discussed in the previous sections. First, it is clear that chi a concentration estimated from the 
blue-to-green band-ratio of ocean color signal in the north polar Atlantic will be significantly 
affected by seasonal shifts in the relationships between absorption and chi a as well as between 
backscattering and chi a. The effects of absorption and backscattering reinforce each other in a 
way that the seasonal difference in the band-ratio chlorophyll algorithm will be quite large. For 
example, Figures 19a and 20a indicate that the relationship for the spring season will predict chi 
a concentration that is higher by a factor of 4-6 compared to that predicted from the summer 
relationship. Second, Figures 19b and 20b suggest that seasonal differences in the band-ratio 
algorithm for estimating total pigments (TChla + AP) will be smaller than those for estimating 
just chi a concentration. Still, estimating total pigments will require a seasonal approach as our 
results suggest that the spring relationship would predict TChla + AP that is higher by a factor 2- 
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4 than the summer estimates. Third, it appears that ignoring seasonal differentiation and using 
just one single algorithm for estimating phytoplankton absorption at 442 nm from the blue-to- 
green reflectance ratio will have a much smaller effect on the accuracy of the resulting product 
compared to pigments (Figures 19c and 20c). The differences in spring and summer estimates of 
a p h(442) are within a factor of 2-2.5. 

Finally, perhaps the most striking and unexpected result is that a single band-ratio algorithm 
based on data pooled together from both spring and summer seasons can provide a fairly good 
estimate of the particulate beam attenuation coefficient at 660 nm (Figures 19d and 20d). Note 
that the relationships of c p (660) versus the product of IOP band-ratio (that is, versus a proxy for 
blue-to-green ratio of reflectance) are the only relationships in Figures 19 and 20, which show no 
clear seasonal shift between the spring and summer data. This result is intriguing and rather 
difficult to interpret in simple terms because the variability in the band-ratio of reflectance in the 
investigated region is driven largely by the variability in absorption and to a smaller extent by 
backscattering [ e.g ., Stramska et ai, 2003] whereas the variability in the beam attenuation in red 
part of the spectrum at 660 nm is expected to be driven primarily by particle scattering, 
especially forward scattering. Although we have already shown that the relatively robust 
relationship between c p (660) and the blue-to-green reflectance ratio in the north polar Atlantic 
can be used as an intermediate step in the regional algorithm for estimating particulate organic 
carbon from ocean color observations [ Stramska and Stramski, 2005], further research is needed 
to examine whether this type of relationship exhibits a similar degree of robustness in other 
regions of the worlds' ocean. 
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Table 1. Input matrix with pigment fractions used in CHEMTAX pigment analysis. 
Abbreviations: Chlorophyll (chi), peridinin (per), 19’butaonoyloxyfucoxanthin (hut), 
fucoxanthin (fuco), 19’hexanoyloxyfucoxanthin (hex), prasinoxanthin (pras), neoxanthin (neo), 
violaxanthin (viola), diadinoxanthin (dda), alloxanthin (alio), lutein, (hit), zeaxanthin (zea). 



chi a 

chib 

chlc3 

per 

but 

fuco 

hex 

pras 

neo 

viol 

dda 

alio 

lut 

zea 

diatoms 

0.541 

0 

0 

0 

0 

0.291 

0 

0 

0 

0 

0.034 

0 

0 

0 

cyano- 

phyceae 

0.615 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0.301 

dino- 

phycaea 

0.308 

0 

0 

0.328 

0 

0 

0.118 

0 

0 

0 

0.074 

0 

0 

0 

prymnesio- 

phyceae 

chryso- 

phyceae 

chloro- 

phyceae 

0.439 

0 

0.13 

0 

0.004 

0.132 

0.187 

0 

0 

0 

0.024 

0 

0 

0 

0.471 

0 

0.057 

0 

0.214 

0.159 

0 

0 

0 

0 

0.027 

0 

0 

0 

0.563 

0.233 

0 

0 

0 

0 

0 

0 

0.032 

0.026 

0 

0 

0.105 

0 

crypto- 

hyceae 

0.75 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0 

0.154 

0 

0 

prasino- 

phyceae 

0.442 

0.363 

0 

0 

0 

0 

0 

0.063 

0.043 

0.014 

0 

0 

0 

0 
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Figure Captions 


Figure 1. Comparison of the total chlorophyll a concentration from HPLC analysis ( TChla ) and 
chlorophyll a concentration estimated with fluorometric and spectrophotometric ( Chla ) methods. 

Figure 2. Average contribution (%) of various algal classes to total chlorophyll a concentration 
{TChla). Polarstem data are shown as averages for the eastern and western sections of the 
transect, as explained in the text. Abbreviations: diatoms ( diat ), prasinophyceae ( prasin ), 
dinophyceae {dino), cyanophyceae {cyan), chlorophyceae ( chi or ), cryptophyceae {crypt), 
prymnesiophyceae ( prymn ), chrysophyceae {chrys). 

Figure 3. Final marker pigment : TChla ratios for selected algal classes, as calculated by 
CHEMTAX. 


Figure 4. Concentrations of (a) photoprotective carotenoids (PPC, see text for definitions), (b) 
photosynthetic carotenoids (PSC), and (c) total accessory pigments (AP) as a function of total 
chlorophyll a concentration {TChla). The results of linear regressions performed on log- 
transformed data are shown as solid lines. The regression line for AP versus TChla from Trees 
et al. [2000b] is shown as the dashed line. 

Figure 5. Approximate contribution (%) of algal size classes to TChla. Abbreviations: micro - 
microphytoplankton; nano - nanophytoplankton, pico - picophytoplankton. 

Figure 6. The absorption coefficient of phytoplankton, a v h, as a function of the total particulate 
absorption coefficient, a v , at 442 and 490 nm. Data from the Oceania and Polarstern cruises are 
indicated by solid and open symbols respectively. 

The least-squares fits for Polarstern data: 
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a ph ( 442) = 0.85 a p ( 442) - 0.002; r 2 = 0.99; n = 73 
a ph (490) = 0.85 a p (490) - 0.001; r 2 - 0.99; n = 73 
The least-squares fits for Oceania data: 
o ph (442) = 0.83 a p (442) - 0.008; r 2 = 0.93; n = 77 
M 49 °) = 0.86 a p (490) - 0.005; r 2 = 0.92; n = 77 

Figure 7. Examples of the spectra of chlorophyll-specific absorption coefficient of 
phytoplankton, o* p h, measured at various chlorophyll a concentrations during Polarstem and 
Oceania cruises. Each example spectrum represents an average obtained from averaging 
absorption spectra from 3 to 5 stations for the range of chlorophyll a concentrations (in mg in" 3 ) 
as indicated on each graph. 

Figure 8. The coefficients of (a) phytoplankton absorption, a p h(A.), (b) total particulate 
absorption, a p (k), plotted as a function of chlorophyll a concentration for the two wavelengths, 
442 and 555 mn. Our data from Oceania and Polarstem are compared with regression lines from 
Bricaud et al. [1995] and Bricaud et al. [1998] shown as solid and dashed lines for 442 and 555 
mn respectively. 

Figure 9. The absorption ratios a w+p h(555)/a w+p h(442) and <7 w +ph(555)/a w +ph(490) as a function of 
TChla from HPLC analysis. Data from Oceania and Polarstem are shown as solid and open 
symbols respectively. 

The least-squares fits for Polarstem data: 
a w+ ph (5 55)/ a w+p h(442) = 1.848 TChla " a541 ; r 2 = 0.93; n = 73 
a w+ph (555)/a w+ph (490) = 1.931 TChla "°- 307 ; r 2 = 0.91; n = 73 
The least-squares fits for Oceania data: 

<2 W + P h(5 5 5)/a w+p h(442) = 1.036 TChla " 0 ' 687 ; r 2 = 0.91; n = 31 
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a w+ph (555)/a w+p h(490) = 1.221 TChla '°- 512 ; r 2 = 0.89; n = 31 


Figure 10. The absorption ratios a w +ph(555)/a w + p h(442) and a w +ph(555)/a w +ph(490) as a function 
of the sum of TChla and AP. Data from Oceania and Polarstem are shown as solid and open 
symbols respectively. 

The least-squares fits for Polarstern data: 

« W + P h(5 55)/ a w+p h(442) = 2.590 (TChla+ AP) ' 0 - 575 ; r 2 = 0.94; n = 73 

a w+ph (555)/a w+ p h (490) = 2.338 (TChla+ AP) ' 0 - 326 ; r 2 = 0.92; n = 73 

The least-squares fits for Oceania data: 

w+ ph ( 555)/ a w+p h(442) = 1.924 ( TChla+ AP) '°- 676 ; ? = 0.93; n = 31 

a w+ p h (555)/a w+ph (490) = 1.945 ( TChla+ AP) ‘°- 512 ; r 2 = 0.93; n = 31 

Figure 11. The absorption ratios a w+p (555)/a w + p (442) and a w +p(555)/<7 w+p (490) as a function of 
chlorophyll a concentration. 

The least-squares fits for Polarstern data: 
a w+p (5 5 5 )/ a w+p (442) = 1.58 TChla ' 0 ' 51 ; r 2 = 0.90; n = 73 
a w+ p(555)/a w+p (490) = 1.77 TChla ' 0 - 3 ; r 2 = 0.88; n = 73 
The least-squares fits for Oceania data: 

<2 w+p ( 5 5 5 )/ a w+p (442) = 0.83 Chla ' 0 - 554 ; r 2 = 0.82; n = 71 
a w+ p(555)/a w+p (490) = 1.04 Chla '°- 427 ; r 2 = 0.79; n = 71 

Figure 12. (a) The total backscattering coefficient bb at 442 and 555 nm versus chlorophyll a 
concentration, (b) Slope parameter y from equation 4 versus the total backscattering coefficient 
bb at 555 nm. 

The least-squares fit for Polarstern data: 
y = 0.026 6 b (555)' 0 ' 750 ; r 2 = 0.95; n = 60 
The least-squares fit for Oceania data: 
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y = 0.067 b h (555)" 0 ' 572 ; r 2 = 0.80; n = 73 

Figure 13. The backscattering ratios 6b(442)/6b(555) and 6b(490)/7>b(555) as a function of 
chlorophyll a concentration. 

The least-squares fits for Polarstern data: 

6b(442)/Z> b (555) = 2.152 TChla - 0093 ; r 2 = 0.70; n = 60 
^ b (490)/6 b (555) = 1.526 TChla - 0052 ; r 2 = 0.73; n = 60 
The least-squares fits for Oceania data: 
b h (442)/b h (555) = 1.583 Chla 4>A06 ; r 2 = 0.66; n = 73 
6 b (490)/6 b (555) = 1.289 Chla - 0 058 ; r 2 = 0.66; n = 73 

Figure 14. The particulate backscattering ratios (a) 6 b p(488)/6 p (488) and (b) 6 b p(660)/6 p (660) as 
a function of chlorophyll a concentration. 

Figure 15. Relationships between the particulate beam attenuation coefficient at 660 nm, 
c p (660), and chlorophyll a concentration. 

The least-squares fit for Polarstern data: 
c p (660) = 0.1995 TChla °' 4995 ; r = 0.89; n = 42 
The least-squares fits for Oceania data: 
c p (660) = 0.4244 Chla °' 5282 ; r 2 = 0.72; n = 60 

Figure 16. The single scattering albedo of particles, co op (488) and co op (660), as a function of 
chlorophyll a concentration. 

Figure 17. The absorption ratios Ow+p(555)/a w +p(442) and a w +p(555)/a w +p(490) as a function of 
the particulate beam attenuation coefficient, c p (660). 

The least-squares fits for Polarstern and Oceania data combined together: 
a w+p (5 55)/ a w+p (442) = 0.6144 c p (660) -°' 6441 ; r 2 = 0.90; n = 102 
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a w+ p(555)/a w+p (490)= 0.3834 c p (660) - 0 - 9065 ; r 2 = 0.89; n = 102 

Figure 18. The backscattering ratios 6b(442)/Z>b(555) and &b(490)/Z?b(555) as a function of the 
particulate beam attenuation coefficient, c p (660). 

The least-squares fits for Polarstern data: 

6 b (442)/6 b (555) = 1.570 c p (660) '° 198 ; r = 0.91; n = 42 

6 b (490)/6 b (555) = 1.283 c p (660) - 0109 ; r = 0.91; n = 42 

The least-squares fits for Oceania data: 
hb(442)/Z> b (555) = 1.426 c p (660) '° 127 ; r 2 = 0.47; n = 60 

6 b (490)/6 b (555) = 1.217 c p (660) - 0 070 ; r 2 = 0.47; n = 60 

Figure 19. Relationships for (a) chlorophyll a concentration versus the product of the 
backscattering ratio and the absorption ratio, R = [6b(442)/6 b (555)] [a w + p (555)/a w + p (442)]. This 
product is an approximation of the blue-to-green ratio of remote-sensing reflectance (see Eq. 1). 
The least-squares fit for Polarstern data: 

TChla = 6.25 17 R ' L5349 ; r = 0.87; n = 42 
The least-squares fit for Oceania data: 

Chla = 1.4437 R ' L32 ; r 2 = 0.83; n = 60 

(b) As in (a) but for the sum of TChla and AP. 

The least-squares fit for Polarstern data: 

(TChla+ AP) = 10.135 R ' L4456 ; r 2 = 0.88; n = 42 
The least-squares fit for Oceania data: 

(TChla + AP) = 4.0145 R - 1 ' 384 ; r 2 = 0.85; n = 26 

(c) As in (a) but for the phytoplankton absorption coefficient at 442 nm, a p h(442). 

The least-squares fit for Polarstern data: 
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a ph ( 442) = 0. 1467 R ' L3449 ; r 2 = 0.98; n = 42 
The least-squares fit for Oceania data: 

<7 ph (442) = 0.0885 R ' 13257 ; r = 0.93; n = 50 

(d) As in (a) but for the particulate beam attenuation coefficient at 660 nm, c p (660). 

The least-squares fits for Polarstern and Oceania data combined together: 
c p (660) = 0.5227 R ' 0 ' 784 ; r 2 = 0.93; n = 92 

Figure 20. Relationships for (a) chlorophyll a concentration versus the product of the 
backscattering ratio and the absorption ratio, R = [6b(490)/6b(555)] [a w+p (555)/a w+p (490)]. This 
product is an approximation of the blue-to-green ratio of remote-sensing reflectance (see Eq. 1). 

The least-squares fit for Polarstern data: 

TChla = 12.3 18 R ' 2 ' 5844 ; r = 0.86; n = 42 
The least-squares fit for Oceania data: 

Chla = 1.6766 R - 17356 ; r = 0.79; n = 60 


(b) As in (a) but for the sum of TChla and AP. 

The least-squares fit for Polarstern data: 

( TChla + AP) = 19.189 R ' 2 ' 4337 ; r 2 = 0.86; n = 42 
The least-squares fit for Oceania data: 

(TChla + AP) = 5. 1763 R - L9446 ; r 2 = 0.85; n = 26 

(c) As in (a) but for the phytoplankton absorption coefficient at 442 nm, a p h(442). 
The least-squares fit for Polarstern data: 

a ph (442) = 0.2696 R ' 2 ' 2766 ; r 2 = 0.97; n = 42 
The least-squares fit for Oceania data: 
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a ph ( 442) = 0. 1054 R - 1 ' 8014 ; r 2 = 0.95; n = 50 

(d) As in (a) but for the particulate beam attenuation coefficient at 660 nm, c p (660). 
The least-squares fits for Polarstern and Oceania data combined together: 
c p (660) = 0.5999 R ' U475 ; r 2 = 0.92; n = 92 
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spectrophotometric or fluorometric Chla [mg m 3 ] fluorometric Chla [mg m‘ 3 ] 


Figure 1. Comparison of the total chlorophyll a concentration from HPLC analysis ( TChla ) and 
chlorophyll a concentration estimated with fluorometric and spectrophotometric (Chla) methods. 
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Figure 2. Average contribution (%) of various algal classes to total chlorophyll a concentration 
( TChla ). Polarstem data are shown as averages for the eastern and western sections of the 
transect, as explained in the text. Abbreviations: diatoms (diat), prasinophyceae (prasin ), 
dinophyceae (dino), cyanophyceae (cyan), chlorophyceae (chlor), cryptophyceae (crypt), 
prymnesiophyceae (prymn), chrysophyceae (chrys). 
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Figure 3. Final marker pigment : TChla ratios for selected algal classes, as calculated by 
CHEMTAX. 
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TChla [mg m' 3 ] 


Polarstem all data 

log(PPC) = 0.816 log (TChla) - 0.948 
r 2 = 0.84, n = 75 

Oceania 

log(PPC) = 0.857 log (TChla) - 0.450 
r 2 = 0.84, n = 31 


Polarstem all data 

log(PSC) = 0.951 log {TChla) - 0.413 
r 2 = 0.96, n = 75 

Oceania 

log(PSC) = 0.947 log (TChla) - 0.085 
r 2 = 0.78, n = 31 


Polarstem all data 

log(AP) = 0.874 log {TChla) - 0.10 

r 2 = 0.98, n = 75 

Oceania 

log(AP) = 0.939 log {TChla) +0.170 
r 2 = 0.87, n = 31 

Trees et al. [2000b] 

log (AP) = 0.934 log (TCh/a)+ 0.028 

r 2 = 0.95, n= 5617 


Figure 4. Concentrations of (a) photoprotective carotenoids (PPC, see text for definitions), (b) 
photosynthetic carotenoids (PSC), and (c) total accessory pigments (AP) as a function of total 
chlorophyll a concentration {TChla). The results of linear regressions performed on log- 
transformed data are shown as solid lines. The regression line for AP versus TChla from Trees 
et al. [2000b] is shown as the dashed line. 
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Figure 5. Approximate contribution (%) of algal size classes to total chlorophyll a concentration 
( TChla ). Abbreviations: micro - microphytoplankton; nano - nanophytoplankton, pico - 
picophytoplankton. 
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Figure 6. The absorption coefficient of phytoplankton, a p h, as a function of the total particulate 
absorption coefficient, a v , at 442 and 490 nm. Data from the Oceania and Polarstem cruises are 
indicated by solid and open symbols respectively. 

The least-squares fits for Polarstern data: 
a ph (442) = 0.85 a p (442) - 0.002; r = 0.99; n = 73 

a ph (490) = 0.85 a p (490) - 0.001; r = 0.99; n = 73 

The least-squares fits for Oceania data: 
a ph (442) = 0.83 a p (442) - 0.008; r = 0.93; n = 77 

M 49 °) = 0.86 a p (490) - 0.005; r 2 = 0.92; n = 77 
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X [ nm ] 

Figure 7. Examples of the spectra of chlorophyll-specific absorption coefficient of 
phytoplankton, a* p h, measured at various chlorophyll a concentrations during Polarstem and 
Oceania cruises. Each example spectrum represents an average obtained from averaging 
absorption spectra from 3 to 5 stations for the range of chlorophyll a concentrations (in mg m" ) 
as indicated on each graph. 
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Figure 8. The coefficients of (a) phytoplankton absorption, a P h(A.), (b) total particulate 
absorption, a p (k), plotted as a function of chlorophyll a concentration for the two wavelengths, 
442 and 555 nm. Our data from Oceania and Polarstem are compared with regression lines from 
Bricaud et a/. [1995] and Bricaud et al. [1998] shown as solid and dashed lines for 442 and 555 
nm respectively. 
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Figure 9. The absorption ratios tf w + P h(555)/a w + p h(442) and a w + p h(555)/a w +ph(490) as a function of 
TChla from HPLC analysis. Data from Oceania and Polarstern are shown as solid and open 
symbols respectively. 

The least-squares fits for Polarstern data: 

«w+ ph (5 55)/ a w+p h(442 ) = 1.848 TChla ‘ 0 ' 541 ; r 2 = 0.93; n = 73 

a w+ph (555)/a w+ph (490) = 1.931 TChla ' 0 - 307 ; r 2 = 0.91; n = 73 

The least-squares fits for Oceania data: 

«w+ ph (5 55)/ a w+p h(442 ) = 1.036 TChla ~ 0M1 ; r 2 = 0.91; n = 31 

a w+ph (555)/a w+ph (490) = 1.221 TChla ' a512 ; r 2 = 0.89; n = 31 
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Figure 10. The absorption ratios a w + P h(555)/a w+p h(442) and a w+p h(555)/a w+p h(490) as a function 
of the sum of TChla and AP. Data from Oceania and Polarstem are shown as solid and open 
symbols respectively. 

The least-squares fits for Polarstern data: 

«w+ ph (5 55)/ a w+p h(442) = 2.590 (: TChla+ AP) '°- 575 ; r 2 = 0.94; n = 73 
a w+ph (555)/a w+ph (490) = 2.338 ( TChla+ AP) '°- 326 ; r 2 = 0.92; n = 73 
The least-squares fits for Oceania data: 

a w+ph (555)/a w+ph (442) = 1.924 (TChla+ AP) ' 0 - 676 ; r 2 = 0.93; n = 31 
a w+ph (555)/a w+ph (490) = 1.945 ( TChla+ AP) '°- 512 ; r 2 = 0.93; n = 31 
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Figure 11. The absorption ratios a w+p (555)/a w+p (442) and a w+p (555)/a w+p (490) as a function of 
chlorophyll a concentration. 

The least-squares fits for Polarstern data: 

n w+p (5 5 5 )/ a w+p (442 ) = 1.58 TChla ' 0 - 51 ; r 2 = 0.90; n = 73 

a w+p (555)/a w+p (490) = 1.77 TChla ' 0 - 3 ; r 2 = 0.88; n = 73 

The least-squares fits for Oceania data: 

n w+p (5 5 5 )/ a w+p (442 ) = 0.83 Chla - 0 ' 554 ; r 2 = 0.82; n = 71 

a w+p (555)/a w+p (490) = 1.04 Chla - 0 ' 427 ; ? = 0.79; n = 71 
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Figure 12. (a) The total backscattering coefficient bb at 442 and 555 rnn versus chlorophyll a 
concentration, (b) Slope parameter y from equation 4 versus the total backscattering coefficient 
bb at 555 mn. 

The least-squares fit for Polarstem data: 
y = 0.026 h b (555)' 0 ' 750 ; r 2 = 0.95; n = 60 
The least-squares fit for Oceania data: 
y = 0.067 b h (555)" 0 ' 572 ; r 2 = 0.80; n = 73 
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0.1 1 10 
TChla or Chla [ mg m ' 3 ] 

Figure 13. The backscattering ratios 6b(442)/Z>b(555) and Z?b(490)/Z?b(555) as a function of 
chlorophyll a concentration. 

The least-squares fits for Polarstern data: 
b h (442)/b h (555) = 2.152 TChla - 0093 ; r = 0.70; n = 60 
6 b (490)/6 b (555) = 1.526 TChla - 0052 ; r 2 = 0.73; n = 60 
The least-squares fits for Oceania data: 

6 b (442)/6 b (555) = 1.583 Chla - 0106 ; r 2 = 0.66; n = 73 
6 b (490)/6 b (555) = 1.289 Chla r 2 = 0.66; n = 73 
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TChla or Chla [ mg m " 3 ] 


Figure 14. The particulate backscattering ratios (a) 6b P (488)/6 p (488) and (b) &b P (660)/& p (660) as 
a function of chlorophyll a concentration. 
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Figure 15. Relationships between the particulate beam attenuation coefficient at 660 nm, 
c p (660), and chlorophyll a concentration. 

The least-squares fit for Polarstem data: 
c p (660) = 0.1995 TChla 04995 ; r = 0.89; n = 42 
The least-squares fits for Oceania data: 
c p (660) = 0.4244 Chla °' 5282 ; r 2 = 0.72; n = 60 
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Figure 16. The single scattering albedo of particles, co op (488) and co op (660), as a function of 
chlorophyll a concentration. 
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Figure 17. The absorption ratios a w +p(555)/a w +p(442) and a w +p(555)/a w +p(490) as a function of 
the particulate beam attenuation coefficient, c p (660). 

The least-squares fits for Polarstern and Oceania data combined together: 

<7 w +p(5 5 5 )/ a w+p (442) = 0.6144 c p (660) -°' 6441 ; r 2 = 0.90; n = 102 
a w +p(555)/a w + p (490)= 0.3834 c p (660) ' 0 - 9065 ; r 2 = 0.89; n = 102 
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Figure 18. The backscattering ratios 6b(442)/Z>b(555) and 6b(490)/Z>b(555) as a function of the 
particulate beam attenuation coefficient, c p (660). 

The least-squares fits for Polarstern data: 

6b(442)/6 b (555) = 1.570 c p (660) ' 0198 ; r 2 = 0.91; n = 42 

6 b (490)/6 b (555) = 1.283 c p (660) -° 109 ; r 2 = 0.91; n = 42 

The least-squares fits for Oceania data: 

6 b (442)/6 b (555) = 1.426 c p (660) ' 0127 ; r 2 = 0.47; n = 60 

6 b (490)/6 b (555) = 1.217 c p (660) - 0 070 ; r 2 = 0.47; n = 60 
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[6 b (442)» 6 (555)] [<.„ p (555)/ a „ p (442)] [6 b (442)» b (555)] [ a ^(555)/a w ^(442)] 


Figure 19. Relationships for (a) chlorophyll a concentration versus the product of the 
backscattering ratio and the absorption ratio, R = [bh(442)//) h (555)] [a w + p (555)/a w + p (442)]. This 
product is an approximation of the blue-to-green ratio of remote-sensing reflectance (see Eq. 1). 
The least-squares fit for Polarstem data: 

TChla = 6.2517 R - 1 ' 5349 ; r 2 = 0.87; n = 42 
The least-squares fit for Oceania data: 

Chla = 1.4437 R ' L32 ; r 2 = 0.83; n = 60 

(b) As in (a) but for the sum of TChla and AP. 

The least-squares fit for Polarstem data: 

(TChla+ AP) = 10.135 R - 1 ' 4456 ; r 2 = 0.88; n = 42 
The least-squares fit for Oceania data: 

( TChla + AP) = 4.0145 R ' U84 ; r 2 = 0.85; n = 26 
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(c) As in (a) but for the phytoplankton absorption coefficient at 442 nm, a p h(442). 
The least-squares fit for Polarstem data: 

a ph (442) = 0. 1467 R ' L3449 ; r = 0.98; n = 42 
The least-squares fit for Oceania data: 
o ph (442) = 0.0885 R ' U257 ; r = 0.93; n = 50 

(d) As in (a) but for the particulate beam attenuation coefficient at 660 nm, c p (660). 
The least-squares fits for Polarstern and Oceania data combined together: 
c p (660) = 0.5227 R '°- 784 ; r = 0.93; n = 92 
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[6 b (490)/6 b (555)][« w+p (555)/a w+p (490)] [6 b (490)/6 b (555)] [a w+p (555)/a w+p (490)] 


Figure 20. Relationships for (a) chlorophyll a concentration versus the product of the 
backscattering ratio and the absorption ratio, R = [6b(490)/6b(555)] [a w +p(555)/a w +p(490)]. This 
product is an approximation of the blue-to-green ratio of remote-sensing reflectance (see Eq. 1). 

The least-squares fit for Polarstem data: 

TChla = 12.3 18 R - 5844 ; r = 0.86; n = 42 
The least-squares fit for Oceania data: 

Chla = 1.6766 R ' 1J356 ; r 2 = 0.79; n = 60 


(b) As in (a) but for the sum of TChla and AP. 
The least-squares fit for Polarstem data: 

(TChla + AP) = 19.189 R ' 2 ' 4337 ; r 2 = 0.86; n = 42 
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The least-squares fit for Oceania data: 

(TChla + AP) = 5.1763 R ' L9446 ; r 2 = 0.85; n = 26 

(c) As in (a) but for the phytoplankton absorption coefficient at 442 mn, n p h(442). 
The least-squares fit for Polarstem data: 

a ph (442) = 0.2696 R ' 2 ' 2766 ; r = 0.97; n = 42 
The least-squares fit for Oceania data: 
a ph (442) = 0. 1054 R ' L8014 ; r = 0.95; n = 50 

(d) As in (a) but for the particulate beam attenuation coefficient at 660 mn, c p (660). 
The least-squares fits for Polarstern and Oceania data combined together: 
c p (660) = 0.5999 R ' U475 ; r = 0.92; n = 92 
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